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Tests of a Habitat Suitability Model For Black-capped Chickadees 


by 


Richard L. Schroeder 


U.S. Fish and Wildlife Service 
National Ecology Research Center 
4512 McMurray Avenue 
Fort Collins, Colorado 80525 


ABSTRACT.—The black-capped chickadee (Parus atricapillus) Habitat Suitability Index (HS1I) 
model provides a quantitative rating of the capability of a habitat to support breeding, based on 
measures related to food and nest site availability. The model assumption that tree canopy volume 
can be predicted from measures of tree height and canopy closure was tested using data from 
foliage volume studies conducted in the riparian cottonwood habitat along the South Platte River in 
Colorado. Least absolute deviations (LAD) regression showed that canopy cover and overstory 
tree height yielded volume predictions significantly lower than volume estimated by more direct 
methods. Revisions to these model relations resulted in improved predictions of foliage volume. 
The relation between the HSI and estimates of black-capped chickadee population densities was 
examined using LAD regression for both the original model and the model with the foliage volume 
revisions. Residuals from these models were compared to residuals from both a zero slope model and 
an ideal model. The fit model for the original HSI differed significantly from the ideal model, whereas 
the fit model for the revised HSI did not differ significantly from the ideal model. However, both the fit 
model for the original HSI and the fit model for the revised HSI did not differ significantly from a 
model with a zero slope. Although further testing of the revised model is needed, its use is recom- 





mended for more realistic estimates of tree canopy volume and habitat suitability. 





The U.S. Fish and Wildlife Service (Service) has pub- 
lished Habitat Suitability Index (HSI) models (Scham- 
berger et al. 1982) for more than 150 species of fish and 
wildlife. These models were developed primarily from 
literature sources and expert reviews and provide a quan- 
titative method to predict habitat suitability based on vege- 
tative and physical characteristics of habitats. Whereas 
earlier Service efforts were directed at model develop- 
ment, the emphasis in recent years has shifted to tests of 
publishe¢ model hypotheses (Schamberger and O’Neil 
1986). I present results from tests of hypotheses in the 
black-capped chickadee (Parus atricapillus) HSI model 
(Schroeder 1983), as well as an empirical test of the HSI 
model itself. 

The black-capped chickadee HSI model provides a 
quantitative rating of the capability of a habitat to support 
breeding, based on measures related to food and nest site 
availability (Fig. 1; Table 1). The food portion of the 
habitat model was developed from data presented in Stur- 
man (1968) and assumes that black-capped chickadee 
population densities are related to arthropod abundance 
and that this is positively correlated with tree foliage vol- 
ume. Sturman (1968) found that canopy volume of trees 
was a strong predictor of black-capped chickadee abun- 


dance (r2 = 0.796). However, measures of tree canopy 
volume as presented in Sturman (1968) are very time 
consuming and may not be suitable for many users of the 
HSI model. An alternative method is presented in the HSI 
model, based on the assumption that measures of tree 
heiGht and tree canopy closure can be used in place of 
Sturman’s volume measurements to assess the suitability 
of the food resource for black-capped chickadees. 

The nest site availability portion of the HSI model as- 
sumes that reproduction, and hence black-capped chick- 
adee numbers, may be limited by the density of snags that 
are 10 to 25 cm diameter at breast height (dbh). This part of 
the model incorporates estimates of the number of suitable 
snags per unit area that are assumed to be required by 
breeding black-capped chickadees. 

The black-capped chickadee HSI model was also tested 
in New York by Bayer and Porter (1988). They compared 
population densities and HSI model values for 28 sites 
using both continuous and discrete (>0.5, >0—0.5, and 0) 
HSI values. The model was successful in predicting popu- 
lation abundance at the discrete level (P < 0.05) but not at 
the continuous level. 

The objectives of this study were (1) to test the hypothe- 
sis that tree canopy volume can be predicted from meas- 
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Fig. 1. Habitat variable and suitability index graphs for the original black-capped chickadee HSI model (Schroeder 1983). 


ures of tree height and canopy closure, and (2) to test the 
relation between the HSI and estimates of black-capped 
chickadee densities during the breeding season. 


Study Area and Methods 


Foliage volume studies were conducted in plains cot- 
tonwood (Populus sargentii) bottomland along the South 
Platte River near Crook, Logan County, Colorado. Black- 
capped chickadee population data were obtained from 
previous studies in the same area conducted by Sedgwick 
and Knopf (1986; J. A. Sedgwick and F. L. Knopf, U.S. 
Fish and Wildlife Service, Fort Collins, Colorado, un- 
published data). Elevations along the South Platte River 
ranged from 1,116 to 1,149 m, and plains cottonwood was 
the dominant overstory tree species. For further descrip- 


tions of the riparian vegetative community refer to Sedg- 
wick and Knopf (1986). 


Volume Test 


A major assumption in the HSI model is that the food 
suitability index (SI) computed through measuring tree 
height and canopy closure is equivalent to the food SI 


Table 1. Life requisite and HSI determination methods 
for the black-capped chickadee (Parus atricapillus) HS/ 
model (Schroeder 1983). 








Model output Determination method 
Food suitability index (SI) = (SIV1 X SIV2)'? or (SIV3) 
Nest SI SIV4 

HSI Lower of food SI or nest SI 














value computed for tree canopy volume as measured using 
Sturman’s (1968) methods. If this assumption of equiva- 
lence is true, canopy volume can be predicted from tree 
height and canopy closure by calculating the food SI for 
tree height and canopy closure and then determining the 
canopy volume associated with the food SI from the graph 
of SIV3 (tree canopy volume). A food SI of 1.0 produced 
through measuring tree height and canopy closure should 
indicate a tree canopy volume 210.2 m3/m2 of ground 
surface. The null hypothesis for the tree canopy volume 
test is that for tree canopy volumes <10.2 m3/m2 ground 
surface there is a linear ratio of 10.2:1 between tree canopy 
volume predicted using the methods of Sturman (1968) 
and the food SI determined through the use of tree height 
and tree canopy closure, with an intercept of 0. 

To test this hypothesis, detailed measurements were 
taken of each tree (minimum height of 5 m) on plots 
selected to provide wide variation in the amount of foliage 
volume estimated to be present. Eighteen plots were se- 
lected with plot sizes of 0.04 ha except for 3 0.1-ha plots in 
scattered mature cottonwood stands and 2 0.01-ha plots in 
dense cottonwood sapling stands. 

Measurements to determine the volume of living fo- 
liage of each tree were taken in May 1988, following the 
methods described in Sturman (1968). Tree height meas- 
urements were taken with a clinometer at a measured 
distance from each tree. The width of individual tree cano- 
pies was measured by vertically projecting the inner and 
outer portion of the living foliage and marking these points 
on the ground with wire flags. 

The percent of tree canopy cover within each plot was 
measured with a spherical densitometer (Lemmon 1956) 
at 20 points, equally spaced along 2 diagonal transects. 
The densitometer was modified to reduce the overhead 
angle viewed in the mirror by using only the four center 
squares. This technique ensured that canopy cover read- 
ings were nearly vertical and did not take in vegetation 
outside the plot. A hit or miss was recorded on each of four 
points within each mirror square. 

Measurements of tree dbh were taken to determine 
basal area for each plot. Basal area measurements were not 
required in the original HSI model. However, basal area is 
widely used and was measured to determine the relation 
between it and tree canopy volume. 


HSI Model Test 


The null hypothesis for the test of the HSI model is that 
regression of the HSI on estimated black-capped chick- 
adee population densities will produce a model with a 
slope significantly different from zero but not signifi- 
cantly different from the slope of a model predicting max- 
imum expected population densities. Black-capped chick- 
adee breeding population data from previous studies 
(Sedgwick and Knopf 1986 and unpublished data) were 
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used to test this hypothesis. These data were collected on 
10 16-ha plots in riparian cottonwood stands in 1982 and 
1984-86, using the variable circular-plot method (Rey- 
nolds et al. 1980). Black-capped chickadee population 
density estimates were determined by using program 
TRANSECT (Burmham et al. 1980). Tree heights were 
measured with a clinometer. Data on snag density, as 
defined in the HSI model, were not collected; however, 
the number of trees with cavities (any opening 3 to 12 cm 
in diameter) was recorded. Based on the scarcity of snags 
on the study site (0.66 snags per hectare), and the fact 
that most cavities (95.8%) were in live trees, it seems that 
the number of trees with cavities provides a reasonable 
substitute for the snag density data required in the HSI 
model. Additional descriptions of methods for population 
estimation and vegetation measurement are presented in 
Sedgwick and Knopf (1986). I obtained measurements 
of tree canopy closure from low-level aerial photographs 
of each of the 10 plots. Individual canopies of cottonwood 
trees were outlined on an acetate overlay of the photo- 
graphs, and digitized. Tree canopy cover was computed by 
comparing the total area of the tree canopies to the area 
of each plot. 


Statistical Analyses 


I used least absolute deviations (LAD) regression (Bar- 
rodale and Roberts 1974) to examine linear relation be- 
tween (1) tree canopy volume predicted from tree height 
and canopy closure, and canopy volume estimated follow- 
ing Sturman (1968), and (2) black-capped chickadee den- 
sities and HSI values. Least absolute deviation regression 
is more resistant to the influence of outlying data values 
than least squares (LS) regression (Wilson 1978; Narula 
and Wellington 1982), and LAD regression estimators 
have smaller variances than LS estimators for data that 
deviate from normality and homogeneous variance, espe- 
cially with small sample sizes (Wilson 1978; Dielman and 
Pfaffenberger 1982). Plots of these data indicated nonnor- 
mality and that heterogeneous variances were present. 
Least squares regression lines and r2 values are given for 
comparative purposes. 

Statistical inferences about proposed and fit regression 
models (Ho: Bfit = proposed; Ha: Bfit + Bproposed) were 
made by using a distribution-free permutation test for 
matched pairs (PTMP; Mielke and Berry 1982) to com- 
pare absolute values of the LAD residuals for alternative 
models. Euclidean distances were used in the PTMP 
analyses because they are consistent with the geometry of 
the LAD regressions (Mielke 1986). The combination of 
LAD regression and distribution-free tests of the LAD 
residuals is a more desirable approach than transforming 
variables to meet LS assumptions because it avoids arbi- 
trarily rescaling data into a geometry that is inconsistent 
with the original measurement scale. 
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Results 


Volume Test 


Least absolute deviation regression lines forced 
through an intercept of 0 (because there can be no volume 
if there are no trees) produced the following fit model 
(Fig. 2): 


Y = 3.21(X), 


where Y = volume (following Sturman’s methods), 
X = food SI, computed as (SIV1 X SIV2)!?. 


According to the HSI model, using Sturman’s (1968) data, 
the proposed model should be of the form Y = 10.2(X). 
Residuals differed between the fit and proposed models 
(P = 0.003), indicating significant differences between 
slopes of the fit and the proposed models. Plots of these 
data (Fig. 2), along with least squares regression analy- 
sis (Y = 3.91[X], r2 = 0.70), indicate that a positive re- 
lation exists between the model output and measured 
foliage volume. 

The plots show that the HSI model assigns high index 
values to plots that have relatively low levels of foliage 
volume. The following changes were made in the model to 
improve the relation (Fig. 3). First, the SI graph for tree 
canopy closure was modified to obtain maximum SI's at 
270% canopy closure. The original model (Fig. 1) pro- 
posed a maximum SI for this variable between 50 and 75% 
canopy closure, with decreasing values above 75%. These 
data indicate that this assumption was not correct for the 
habitat type studied. The modified relation seems to more 
accurately depict the effects of changes in canopy closure 
on foliage volume. 
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Fig. 2. Plots of original food suitability indices for cottonwood 
habitats versus tree canopy volume measured following Stur- 
man (1968). . 


Second, the original variable “‘average height of over- 
story trees’’ considered only the height of overstory trees 
(defined as trees 280% of the height of the tallest tree, 
where trees are considered to be any woody vegetation 25 
m in height). Thus, the indirect estimates of canopy vol- 
ume in stands with a few very tall trees and many trees 
<80% of the height of the tallest tree would be overrated 
with the original modei because the limited volumes of the 
many smaller trees would not be averaged with the taller 
trees. Therefore, the model was revised to consider the 
average height of all trees. The original model attained 
maximum SI values at tree heights >15 m. These data 
indicate that volume continues to increase as tree heights 
increase above 15 m. Therefore, the SI relation was mod- 
ified to provide increasing values as tree heights increased 
up to 25 m. In addition, the original model provided SI 
values for woody vegetation of any height >0 m. Thus, a 
stand with an average tree height of 5 m would receive an 
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SI of 0.33. These stands would not provide suitable breed- 
ing habitat for the black-capped chickadee, which is a 
forest-dwelling bird. Therefore, the SI relation was mod- 
ified to provide SI values of 0 for stands with average tree 
heights of $5 m. 

The last change made in the model was to eliminate the 
square root function from the food SI calculation. The 
volume of a three-dimensional object is equal to the prod- 
uct of the area of the base and the height. For a forest 
stand, a measure of tree canopy closure is comparable to 
the area occupied by tree canopies, and a measure of 
average tree height considers the vertical height. Thus, it is 
more logical to compute the food SI by using the product 
of SIV1 (canopy area) and SIV2 (height), rather than the 
geometric mean. 

These model revisions had a marked effect on the abil- 
ity of the model to predict tree foliage volume. Least 
absolute deviation regression analysis of the revised food 
SI relations, forcing a 0 intercept, produced the following 
fit model (Fig. 4): 


Y = 10.11(X), 


where Y = volume (following Sturman’s methods), 
X = food SI, computed as (Rev. SIVI X Rev. 
SIV2). 


Residuals from the fit and proposed (Y = 10.2[X]) models 
did not differ (P = 0.61), indicating no significant differ- 
ences between slopes of the fit and proposed models. The 
least squares regression fit, again with intercept forced to 
0, was Y = 10.39(X) (r2 = 0.877). 

Data on basal area for each plot were also regressed 
against volume as computed using Sturman’s (1968) 
methods, and it seems that basal area is a good predictor of 
foliage volume. Least absolute deviation regression analy- 
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Fig. 4. Plot of revised food suitability indices for cottonwood 
habitats versus tree canopy volume measured following Stur- 
man (1968). 
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sis, with the intercept forced to 0, produced the following 
fit model: 


Y = 0.138(X), 


where Y = volume (following Sturman’s methods), 

X = basal area from plot data (m2/ha). 
Residuals from the fit models and models with a slope of 
zero were different (P = 0.05), indicating the slope of the 
fit model was significantly different than zero. 


HSI Mode! Test 


The test of the volume relation noted previously only 
considers the model’s ability to predict tree canopy vol- 
ume based on measures of tree height and canopy closure 
or basal area. A significant additional question is the rela- 
tion of tree canopy volume to measures of black-capped 
chickadee abundance. The hypothesis in the HSI model is 
that habitats will provide increasing potential black- 
capped chickadee population densitics as the HSI in- 
creases from 0 to 1. 

Both the original and the revised food SI relations were 
used to test the relation between the HSI and black-capped 
chickadee population densities from the 10 16-ha riparian 
cottonwood plots. The model states that the HSI equals the 
lower of either the food or nest SI values. Population 
densities, original and revised food SI's, nest SI’s, and 
original and revised HSI values for the 10 plots are pro- 
vided in Table 2. Least absolute deviation regressions of 
population densities on the original and revised HSI val- 
ues produced the following fit models: 


Original HSI Model: 
Y = 0.098 + 0.416(X), 


where Y = population density of black-capped chick- 


adees (both sexes) per hectare, 
X = original HSI. 
Revised HSI Model: 


Y = 0.096 + 2.344(X), 


where Y = population density of black-capped chick- 
adees (both sexes) per hectare, 
X = revised HSI. 


Residuals of the fit models and a model with a zero slope 
did not differ (original HSI model, P = 0.20; revised HSI 
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Table 2. Black-capped chickadee (Parus atricapillus) densities, food and nest SI values, and HSI values for the original and revised 








models. 
Plot Chickadees Original Revised Nest Original Revised 
number per hectare food SI food SI SI HSI HSI 
l 0.23 0.62 0.097 0.88 0.62 0.097 
2 0.39 0.70 0.084 0.81 0.70 0.084 
3 0.46 0.62 0.080 9.63 0.62 0.080 
4 0.37 0.81 0.102 0.57 0.57 0.102 
5 0.30 0.58 0.087 0.39 0.39 0.087 
6 0.30 0.68 0.113 0.41 0.41 0.113 
7 0.29 0.60 0.071 0.58 0.58 0.071 
8 0.10 0.58 0.048 0.19 0.19 0.048 
9 0.14 0.47 0.052 0.10 0.10 0.052 
10 0.15 0.49 0.021 0.31 0.31 0.021 





model, P = 0.34), indicating that neither of the fit models 
differed significantly from a zero slope. 

An important aspect of testing HSI models is to com- 
pare the model output to an estimate of maximum ex- 
pected performance of the species. The black-capped 
chickadee HSI model was developed to encompass the 


its maximum. The LS regression fit for the revised HSI 
model was Y = 0.063 + 2.772(X) (r2 = 0.425), yielding a 
maximum estimated population density of 2.84/ha when 
the HSI is 1. 


entire range of the species. To determine the highest popu- Discussion 
lation densities that have been encountered in previous 
studies, data from more than 350 Breeding Bird Census 

Volume Test 


(BBC) reports (T. Engstrom, Cornell Laboratory of Or- 
nithology, Ithaca, New York, personal communication) 
were analyzed. These data indicated that the average den- 
sity of the highest five BBC’s for black-capped chickadees 
was 1.4 males per hectare (excluding one small orchard 
with | breeding male and a corresponding density of 2.5 
males per hectare). Assuming that actual density of both 
sexes would be double this figure (2.8 chickadees per 
hectare), a proposed model for an ideal HSI would have 
the form 


Y = 2.8(X), 


where Y = population density of black-capped chick- 
adees (both sexes) per hectare, 
X = HSI. 


Residuals from the fit model for the original HSI and the 
proposed ideal model were different (P = 0.002), indicat- 
ing a significant difference between the slopes of these 
models. This difference can be seen because the fit model 
for the original HSI yields a predicted population density 
of 0.52/ha when the HSI is at its maximum value of |. Re- 
siduals from the fit model for the revised HSI and the pro- 
posed ideal model did not differ (P = 0.50), indicating no 
significant difference between the slopes of these models. 
The fit model for the revised HS! predicts a population 
density of 2.44/ha when the HSI is at 1. Thus, the revised 
HSI model more closely approximates maximum ex- 
pected population densities (2.8/ha) when the HSI is at 


Proper evaluation of HSI models involves assessing the 
performance of both the HSI and other testable assump- 
tions within the model. Schamberger and O’ Neil (1986) 
noted that most studies tested the HSI against some meas- 
ure of population status (Clark and Lewis 1983; Cook and 
Irwin 1985; Lancia and Adams 1985). Here, I have as- 
sessed the performance of both a model assumption and 
the HSI value. 

The foliage volume portion of the tests indicates that the 
original assumed relation of canopy cover and overstory 
tree height to foliage volume, although showing a strong 
linear relation, yielded predictions significantly lower 
than the volume estimated by Sturman’s (1968) methods. 
Revisions to these model relations resulted in improved 
predictions of foliage volume. 

Habitat Suitability Index model testing can lead to new 
information if the test is not limited to hypotheses con- 
cerning only the variables found in the HSI model. The 
measurement of tree dbh was added to this study to assess 
the possibility that basal area could be used to predict tree 
foliage volume. Basal area data are widely available and 
are often used in models to predict future forest condi- 
tions. Thus, it would be a useful variable for biologists 
conducting forest habitat studies. Data for the cottonwood 
habitats indicate that basal area is an excellent predictor of 
foliage volume. Based on these results, users of the black- 
capped chickadee HSI model have three reasonable op- 
tions to estimate foliage volume: Sturman’s (1968) 














methods, the revised SI graphs and combination formula 
for tree height and canopy closure, and the predictive 
model using basal areu. 

Sturman (1968) provided only two potential tree saapes 
to be used in calculating foliage volume for an individual 
tree, one each for conifers and deciduous trees. Many 
cottonw dod trees do not clearly match these shapes. Given 
the irregular shapes of many trees, | made visual estimates 
and adjustments of the outer and inner ~anopy layers to 
approxiinate the desired shape. Mawson et al. (1976) de- 
scribed a method to determine tree crown volume using 15 
tree shapes. Improved estimates of volume may be possi- 
ble by expanding the number of tree shapes to match more 
realistically what is encountered in forested habitats. 


HSI Model Test 


Results of the test of the HSI as a predictor of black- 
capped chickadee population density are not as easily 
interpreted as the volume test results. The fit model for the 
revised HSI did not differ from the ideal model, whereas 
the fit model for the original HSI differed significantly 
from the ideal model. However, both the fit model for the 
original HSI and the fit model for the revised HSI did not 
differ significantly from a model with a zero slope. Given 
the changes made in the revised HSI model and the spe- 
cific data set used to test the model, several possible expla- 
nations can be offered for these results: (1) Whereas the 
revised model seems to be a better predictor of tree canopy 
volume, volume may be an inadequate predictor of popu- 
lation densities in cottonwood sites. Additional factors 
may need to be added to the model to improve its predic- 
tive ability; (2) the number of sites (n = 10) was inadequate 
to effectively test the HSI model. If the same spread as in 
these data occurred over 30 sites, the revised fit model 
would differ significantly from a zero slope (P = 0.047), 
and thus fail to reject both portions of the null hypothesis. 
In this same situation (if n = 30) the original fit model 
would continue to be significantly different from the pro- 
posed ideal model (P < 0.0001), thus rejecting this portion 
of the null hypothesis; and (3) the 10 sites used for the 
black-capped chickadee test had too narrow a range of tree 
canopy cover values (11 to 33%) and average height of 
trees (7.5 to 12.3 m) to provide an adequate test of the 
hypothesis. Without further study, it is not possible to 
determine which of these possibilities best explains the 
results obtained in this test. 

Only a single measure of black-capped chickadee per- 
formance (number per hectare) was used in this test. How- 
ever, population density may be a misleading indicator of 
habitat quality, as noted by Van Horne (1983). Additional 
measures of performance, such as nest success or number 
of young fledged, might improve the test. In addition, the 
HSI, as a measure of key habitat variables, considers only 
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a subset of all of the factors that may influence individuals 
and resultant populations. Thus, habitat measurements 
may indicate ideal conditions, but nonhabitat-related fac- 
tors (such as disease or severe weather) may cause popula- 
tion declines. The HSI is best viewed as an index of 
habitat-imposed limitations on a population. For example, 
an HSI of | would indicate no habitat-related limitations 
exist, whereas decreasing HSI values would indicate in- 
creasingly significant habitat limitations. Actual popula- 
tion densities should not exceed the levels indicated by the 
HSI values but may range below the HSI predictions and 
not falsify this hypothesis. Specific tests related to this 
aspect of the HSI have been suggested (U.S. Fish and 
Wildlife Service 1987, 1989), but the sample size in this 
study was too small to conduct the tests. 

One suggested change in the nest portion of the black- 
capped chickadee HSI model resulted from this study. The 
original model required a measure of the number of snags 
from 10 to 25 cm dbh as an indication of nest site avail- 
ability. Sedgwick and Knopf (1986), however, noted that 
snag densities alone do not accurately reflect the availabil- 
ity of suitable nesting substrates in riparian cottonwood 
stands. The majority of potential nest sites in these stands 
occurred in cavities in live trees. In eastern hardwood— 
coniferous forests, use of live trees is much less common 
(of 44 nests, 41 were in snags or stubs and only 3 were in 
live trees; Runde and Capen 1987). The best overall meas- 
ure of black-capped chickadee nest site availability may 
be the combined density of the number of trees with 21 
cavity (minimum diameter of 10 cm) and the number of 
snags (minimum dbh of 10 cm). Further studies are re- 
quired to test whether regional differences in use of live 
trees, as opposed to snags, should be incorporated into the 
HSI model. 

I have not developed alternatives to the revised black- 
capped chickadee HSI model, nor have I attempted to 
search out the best model (i.¢., best statistical fit) from the 
data. Fitting anew model to these 10 data points may result 
in a product with no better predictive ability than the 
current revision. I recommend that additional tests be per- 
formed on the revised HSI model before further changes 
are incorporated. In the interim, I recommend that the 
revised HSI model be used by field biologists and wildlife 
managers. The revised model uses the SI relations for 
canopy closure and tree height illustrated in Fig. 3, the 
food SI computed as the product of SIV1 and SIV2, and 
the variable for nest suitability (SIV4) measured as the 
combined density of the number of trees with 21 cavity 
(minimum diameter of 10 cm) and the number of snags 
(minimum dbh of 10 cm). Although additional testing is 
necessary, it seems that the revised model will provide 
more realistic estimates of both tree canopy volume and 
black-capped chickadee habitat suitability than the origi- 
nal version. 
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